Production efficiency of Feshbach molecules in fermion systems 
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We present a consistent nonequilibrium theory for the production of molecular dimers from a 
two-component quantum-degenerate atomic Fermi gas, via a linear downward sweep of a magnetic 
field across a Feshbach resonance. This problem raises interest because it is presently unclear as 
to why deviations from the universal Landau-Zener formula for the transition probability at two- 
level crossing are observed in the experimentally measured production efficiencies. We show that 
the molecular conversion efficiency is represented by a power series in terms of a dimensionless 
parameter which, in the zero-temperature limit, depends solely on the initial gas density and the 
Landau-Zener parameter. Our result reveals a hindrance of the canonical Landau-Zener transition 
probability due to many-body effects, and presents an explanation for the experimentally observed 
deviations [K.E. Strecker et al, Phys. Rev. Lett. 91, 080406 (2003)]. 
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The advances of last years in the experimental tech- 
niques of atomic and molecular trapping and coohng, 
combined with the possibility of externally tuning the 
inter-atomic interactions, have ushered in a series of novel 
applications: the emergence of a molecular Bose-Einstein 
condensate (BEC) from a Fermi gas [^Q, observation of 
coherent oscillations between an atomic condensate and 
molecules , formation and propagation of matter- wave 
soliton trains , and the examination of Cooper pairing 
in the BCS-BEC crossover regime 

Regardless of their role in further investigations, in 
all these experiments diatomic molecules are being pro- 
duced from an atomic BEC or a quantum degenerate 
two-component Fermi gas either by sweeping a magnetic 
field across a Feshbach resonance (FR) IslTol polllTLIl^ or 
by two-photon Raman photo-association |l2l |. Therefore, 
the understanding and control of the molecular produc- 
tion mechanism is of special importance. 

Notwithstanding the differences in the details of FR 
experiments, they all show a growth of the molecular 
conversion efficiency (MCE) with the inverse sweeping 
rate of the magnetic field, , that saturates at values 
less than 100% in the adiabatic regime. 

The attempts aimed at explaining the dependence 
of MCE on B, resonance width, initial atomic den- 
sity and temperature for two-component Fermi sys- 
tems can be broadly classified into two classes: i)semi- 
phenomenological scenarios [T3. ITEI that reduce the 
many-body physics to a two-atom description modeled as 
a two-state Landau-Zener (LZ) system [l| corresponding, 
respectively, to the free two-atom scattering state and the 
bound molecular state, and ii) numerical many-body cal- 
culations 0,^3 based on an effective Hamiltonian first 
proposed by E. Timmermans et al. |l7j. 

Class i) is appealing by its use of simple and intu- 
itive physical pictures, but their predicted (temperature 
independent) upper MCE limit of 50% contradicts the 
experimentally observed far greater values 0, ll^. The 
recent experimental work by E. Hodby et al. |l3l | also 



shows a pronounced T-dependence of this upper limit. 
The breakdown of the simple two-level LZ picture can 
be corrected only by introducing suplimentary ad hoc 
assumptions in these semi-phenomenological scenarios, 
whereas it emerges naturally from a bona fide many-body 
analysis (see below). The work in Class ii) has shown, 
albeit under some simplifying assumptions, the potential 
of Hamiltonian 17] in analyzing the temperature depen- 
dence of the MCE saturation in the adiabatic regime. 

In this Letter we study the atom-molecule conversion 
in ultra-degenerate two-component Fermi gases subject 
to a linear downward sweep of a magnetic field across an 
s-wave FR, in the spirit of experiments |^ ITol ITll [l3| . 
We focus on the zero-temperature dynamics, and de- 
velop a nonequilibrium theory, pertinent to both weak 
and strong atom-molecule coupling (measured in Fermi 
energy units), which allows for a full account of the ef- 
fects of quantum statistics. The MCE is calculated in 
terms of real-time Green functions (GF), and represented 
as a power series in terms of a dimensionless parame- 
ter that depends only on the initial gas density and the 
LZ parameter. An exact evaluation of Feynman-Keldysh 
diagrams for second and fourth order processes reveals 
a clear deviation from the LZ transition probability at 
two-level crossing. This deviation, whose origins reside 
solely in many-body effects, signals a suppression of the 
LZ-predicted MCE even for moderately small values of 
-B~^, as observed experimentally in |lOj] . Equally impor- 
tant, our MCE result doesn't display an a priori upper 
limit of 50% at T = as suggested in |0], though fur- 
ther work is necessary to establish the correct limit and 
its T-dependence [Tsj] . 

The starting point of our analysis is the Hamiltonian 
[itL describing a system of fermionic atoms FR- 
coupled to bosonic molecules, H{t) ~ Ho{t) + V, with 
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where a^(p), a{p) and w{p), h{p) are fermionic creation 
and annihilation operators describing atoms of momen- 
tum p and "spins" t (a) and J, (5), respectively, and P{p), 
f{p) play the same role for the bosonic molecules. Other 
quantities entering H are e^{p, i) — e^,(p)~ ^^B(t), with 

— a,b, /, where /i^ is the projection of the magnetic 
moment along the direction of the magnetic field B{t) 
with which interacts via Zeeman coupling, and e^{p) 
is the dispersion relation which accounts for the single- 
particle energy renormalization due to nonresonant col- 
lisions, and simply reduces to the kinetic energy p^/2m^ 
in a coUisionless regime g is the two- atom-molecule 
coupling [23| which controls the FR width and V is the 
volume of system. The free two-atom scattering state and 
the molecular state (MS) have different spin configura- 
tions and their coupling is mediated via the intra-atomic 
hyperfine interaction [21| which flips the electronic and 
nuclear spins of one of the colliding atoms. Depending 
on the magnetically tuned energy difference between the 
two states, the MS is quasi-bound (virtual) and belongs 
to a closed scattering channel if its energy exceeds that 
of the two-atom channel, becomes resonant with the lat- 
ter when their energies are equal, and turns truly bound 
when its energy is the lesser of the two. 

In order to probe the MCE dependence on B we eval- 
uate real-time GF within the Keldysh formalism (KF) 
psf . The method is based on the use of a closed contour 
for time ordering, which runs from — oo to +oo and then 
back to —oo. Both branches of the contour propagate 
along the real time axis and any point along them can 
be characterized by two parameters, written compactly 
as t'^ , with T being the time variable and 7 a bookkeep- 
ing index that distinguishes between the forward (7 = -I-) 
and reverse (7 = — ) time directions. The basic quantities 
of KF are the contour-ordered real-time GF: 
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with Q — A, B, T ior 'ijj = a, b, /, respectively, a, 
/3 = ± and ((• ■ ■)) = Tr[p(to) (•••)]; Pi^o) is the initial 
density operator at — —00, are the Heisenberg- 
picture (HP) operators relative to ioj and Tc is a contour- 
ordering operator. The corresponding free GF read 
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iQo {Pl,Tl]P2,T2) = 

where ipi are the interaction-picture (IP) operators rela- 
tive to tQ. 

In experiments 0, 0, 0, 0| an ultracold two- 
component Fermi gas is prepared as an incoherent mix- 
ture of equal populations in each state, and extreme 
quantum-degeneracy, at temperatures as low as T ^ 
O.ObTp, has been reached where Tp is the Fermi 
temperature. In this regime the fall-off of the Fermi dis- 
tribution from 1 to takes place in an extremely narrow 



energy interval '--^ O.OSef, where ep is the Fermi energy, 
and the fuzziness of the Fermi surface becomes virtually 
unimportant. In this vein we take p{to) = |$o) (^o|) 

with |$o) = ^a^p)b^p)\VAC), where |VAC) is the vac- 
uum state [25j . The average number of molecules at time 
t is given by 



(4) 



and, upon expressing the tpn operators in terms of their 
IP form ?/>/, ipH{t) = u}{t,to)ipi {t)Ui{t,to), and ex- 
panding the IP time-evolution operator Ui{t, to) as a for- 
mal series in the coupling constant g, the following sys- 
tematic expansion of iJ^~^ (je,t;k,t^ ensues: 



/+00 /•+00 
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where Vi is the IP form of V , and the sum X]{7}=± runs 
over all n-tuples (71, . . . , 7„) with 7^ = ±. 

Since // (p, t) |$o) ~ for any p and t, and Vj ^ fi + 
/|, it follows that only terms with even n can have a 
nonvanishing contribution in Eq. 

Due to the form of |<i>o), the free GF 
Gq'^ {pi,Ti;p2,T2) cc S{pi,p2), where S{pi,p2) is 
the Kronecker delta, and their expressions ,26:j are 

zg+- {p- n , r2) - -6 {ep ~ {p}) (6) 
(p; n, r2) - e {e.^{p) - ep) " -^ip-^)^^^ (7) 



with Q — A, B for ip = a, 6, respectively, and 
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and finally (p; n, T2) = 9 (x) 5^+ {p;ti,T2) + 

e{-x)g+- {p;ti,T2), g (p;ti,T2) = 

0{x)g+^ {■p;Ti,T2) + 9 {-x)g-+ {■p;Ti,T2) for any 
g = A, B, where x — ti — T2 and 9 (x) is the 
Heaviside function. 

Each average corresponding to the terms in Eq. can 
be performed by means of a generalized version of Wick's 
theorem in which the contractions are defined with re- 
spect to the contour-ordering operator Tc, and a dia- 
gram is associated with each way of contracting the field 
operators into pairs These diagrams have the same 
topology as those occurring in the ordinary quantum field 
theory (OQFT) for systems in equilibrium |27j . the only 
difference being an additional label 7 = ± that has to 
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the contribution from the (D2) diagrams is 
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FIG. 1: Feynman-Keldysh diagrams for second (Dl) and 
fourth order (D2 - D4) processes. The free Green functions 
are represented by continuous Unes for a-fermions (^Iq''''), by 
dashed lines for 6-fermions (Z3q'^^), and by wiggly lines for 
bosons (jrj'''^ ). 



be attached to each interaction vertex. As in OQFT, the 
disconnected diagrams corresponding to vacuum polar- 
ization vanish ,2;||, and only topologically distinct dia- 
grams need to be considered. 

Since we are interested in analyzing the dependence 
of MCE on , and not the behavior of the average 
number of molecules in time, we set the initial time of 
atomic gas preparation at = —00, and the molecule- 
counting time at tm — 00. The diagrams representing 
the contribution from second and fourth order processes 

to l^^f'^ shown in Fig. The (Dl) diagrams 

contribute as 
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and each of the diagrams (D3) and (D4) contributes 
equally as 
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pi,qi P2 
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dTi - ■ ■ / dT4 iT^^ (pi -f (fi ; 00, Ti) 
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where Ml = (/iy — fia ^ ^-bjB^ Nq is the total number of 
atoms present in the system before the magnetic field is 
applied, and hq = Nq/V is the initial density. 

The evaluation of all integrals entering Eqs. Hl()|l — 112(1 
can be carried out exactly, without any supplementary 
assumptions, and a presentation of the steps involved is 
beyond the scope of this letter. 

Upon collecting results and introducing the notation 

F = 27raz(VJ|), where 
LZ parameter [l|, we obtain 



is the canonical 



MCE = 



2(iV/)(oo) 
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The n-th term of this series is represented by the set 
of Feynman-Keldysh diagrams containing 2n vertices. A 
similar approach, i.e. the development of a formal solu- 
tion in terms of an infinite scries in powers of the cou- 
pling constant, with further calculation of every term and 
summation of the resulting algebraic series, has been em- 
ployed to exactly calculate the probability of nonadia- 
batic transitions in a multiple-crossing LZ model and 
it constitutes a powerful alternative to Landau's method 
Pi of analytic continuation in complex time. 
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Eq. (|13|l reveals deviations from the universal two- 
level LZ formula and also from the phenomenological 
correction proposed in as 

v{i-e~^)=v(r-lr' + o{r')y (i4) 

where r] < 50% is a constant that depends on the initial 
populations of the two-component Fermi gas. 

Since > ^, Eq. (fT!^ shows that, as increases, 
the MCE grows slower then predicted by the LZ formula, 
and this behavior is experimentally supported [lol |. In 
our theory, the approach towards saturation is not due to 
a mere contraction of the LZ formula by a multiplicative 
factor determined solely by the initial state preparation, 
as proposed in the LZ scenarios 0, , but has a rather 
dynamical nature as the atom-molecule conversion takes 
place in a many-body medium in which the effects of 
quantum statistics play a crucial role. 

Examination of higher order diagrams indicates that 
MCE is a function depending solely on the parameter F. 
Therefore, in the extreme adiabatic regime, correspond- 
ing to F — > oo, MCE must have a universal limit at T = 
which, unlike in the phenomenological result (I14|l . is not a 
priori bounded by 50%. In practice, as the experiments 
are carried out at finite T and F, the smearing of the 



Fermi surface when T approaches Tp, and the quantum 
degeneracy reaches its lower limit, must be taken into 
account |24| for analyzing the T-dependence of the MCE 
saturation 'l3l|. 

In conclusion, we have considered the MCE for a 
hyperfine-induced s-wave FR and developed a consistent 
many-body nonequilibrium theory, based on the real- 
time GF approach, in which all atomic and molecular 
states are included, and the effects of quantum statistics 
are fully accounted for. 

We demonstrated, by analytically evaluating the MCE 
up to fourth order in the hyperfine coupling constant, 
that the canonical LZ formula at two-level crossing is vi- 
olated in this system due to many-body effects which sys- 
tematically decrease the LZ transition probability, even 
for moderately small values of B~^. This result indicates 
that in degenerate Fermi gases the effects of quantum 
statistics near a FR play a crucial role, and the singling 
out of independent two-atom pairs from an ensemble of 
delocalized indistinguishable particles, as proposed in the 
LZ scenarios untenable. 
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